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Abstract. For a class of dynamical systems, the cookie-cutter maps, we prove 
that the Lyapunov spectrum coincides with the map given by the Newton- 
Raphson method applied to the derivative of the pressure function. 



1. Introdcution 

In this note we establish a relation between the Newton-Raphson method for 
finding roots of an equation and an important function in the dimension theory of 
dynamical systems, namely the Lyapunov spectrum. 

The dimension theory of dynamical systems has received a great deal of attention 
over the last years (see [Bal IPej ) . Multifractal analysis is a sub-area of dimension 
theory devoted to study the complexity of level sets determined by invariant local 
quantities. Usually, the geometry of the level sets is complicated and in order to 
quantify its size the Hausdorff dimension is used. In this note we will focus our 
attention in one of these local quantities, namely the Lyapunov exponent. Let 
T : I — > I be a piecewise differentiable map defined on the unit interval I = [0, 1]. 
The Lyapunov exponent of T at the point a; € J, is defined by 

X T (x) = X(x) = lim -log\(T n Y(x)\, 

n— >oo 71 

whenever the limit exists. It is a dynamical quantity that measures the exponential 
rate of divergence of infinitesimally close orbits. It is possible for the Lyapunov 
exponent to attain a wide range of values. In order to understand the complexity 
of the induced level sets we study the Lyapunov spectrum, which is the function 
defined by 

L(a) := dimjj({a; G / : X(x) = a}). 

Here, dim# denotes the Hausdorff dimension of a set (see |Fal Chapter 2] for a 
precise definition). This function can be described using thermodynamical for- 
malism, in particular it can be written as a function of the topological pressure 
P{— t log |T"|), see Section [3] for precise definitions and statements. 

We will relate the Lyapunov spectrum with the function used to apply the 
Newton-Raphson method. Given a real function / : R — > M, under certain assump- 
tions, the iterates of the Newton map Nf(x) converge to a root of the equation 
f(x) — (for precise statements see section [2T2|) . The study of this function and 
its dynamic behaviour has attracted a great deal of attention (see [MiJ. 
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The main result in this note is that, for a certain class of dynamical systems, the 
Lyapunov spectrum equals the Newton map of the topological pressure P(—t log \T'\), 
that is 

L(-P'(-t log \r\))=N P (t). 



2. The Legendre transform and the Newton method 

Let / : M — > M be a strictly decreasing convex function of class C 2 . In what 
follows we will dchnc its corresponding Legendre transform and the associated 
Newton map. Moreover, we will show how are they related. 

2.1. Legendre transform. Let t' G R, the line S a = a(t — t') + f(t') is called 
a support line at t' if it always lie below the graph of /. There exists an interval 
of the form [a m ; n , a max ] (this interval need not to be bounded) such that if a G 
[dmini <*max] the graph of / has a support line S a of slope equal to —a and for 
such a the support line is unique. The Legendre transform of / is the function 
F : [a m \ n , a max ] — > R defined as the value of the intercept of S a with the vertical 
axis. That is 

(1) F(a) = m£{f(t) + at : t G R}. 

It turns out that if t a G R is the unique point for which f'(t a ) = —a then 

(2) F{a) =f(t a ) + at a . 

Let us remark that in this simple setting the support line S a is the tangent to the 
graph of / at the point (t a ,f(t a )) and that the interval [a m i n ,a max ] is uniquely 
determined by the range of the derivative of /. 

Geometrically, the Legendre transform is the intercept of the tangent line of / 
with the vertical axis. 



2.2. The Newton-Raphson method. The Newton-Raphson method is an al- 
gorithmic method that provides a sequences of approximations to the root of the 
equation f(t) = 0. The algorithm is based on the study of the dynamics of the so 
called Newton map. This map is defined by 

(3) N f (t)=t-M. 

Note that since the function / is strictly decreasing the Newton map is defined over 
the real numbers (the function / has no critical values). It turns out that a point 
d G R is a root of / if and only if it is a fixed point for Nf. Moreover, since 

(4) - {m) 2 . 

we have that N'(d) = 0. In particular, the point d is an attractor for Nf and 
iterations of the Newton map converge to the root d. That is, for any x G R the 
sequence {x, Nf(x), Nj(x), . . . , Nj(x), . . . } converges to the root. 

Geometrically, the procedure involved in this algorithm is to consider the tangent 
to the graph of / at the point (t, f(t)) and define the Newton transformation as 
the intercept of the tangent with the horizontal axis. 
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2.3. A relation between the Legendre transform and the Newton method. 

In virtue of the above we have that in order to compute the Legendre transform 
and the Newton map of / we need to compute the tangent to /. 

Lemma 2.1. Let f : R — > R be a strictly decreasing convex function of class C 2 
and let t a el be a point such that f'(t a ) = —a. Then 

(5) N f (t a ) = ±F(a). 

Proof. Note that since t a G R is a point such that f'(t a ) = —a then in virtue of 
equation ([3]) the Newton map is given by 

N f (t a ) = t a +^A 
a 

and the Legendre transform (see equation @) is given by 

F(a) = f(t a ) + at a . 

The result now follows. □ 

3. Thermodynamic formalism 

A large class of interesting dynamical systems have many invariant measures. It 
is, therefore, an important problem to find criteria to choose relevant ones. Ther- 
modynamic formalism is a set of ideas and techniques which derive from statistical 
mechanics and that was brought into dynamics in the early seventies by Ruelle and 
Sinai among others (see the books [Kj IRu2[ IWaj for a detailed and nice exposition 
of the subject). This formalism provides procedures for the choice of interesting 
measures. In this note we will consider the following class of dynamical systems. 
Given a pairwise disjoint finite family of closed intervals Ji,...,7 n contained in 
[0, 1] we say that a map: 

n 

T : |J h -> [0, 1] 

»=i 

is a cookie-cutter map with n branches if the following holds: 

(1) T{h) = [0, 1] for every i £ {1, . . . , n}, 

(2) The map T is of class C 1+£ for some e > 0, 

(3) \T'(x)\ > 1 for every x G h U • • • U /„. 

The dynamics is concentrated in the repeller A of T, which defined by 

oo 

A := p) T~ n I. 

71 = 

One of the most important objects in thermodynamic formalism is the so called 
topological pressure of T. It is defined by 

(6) P{t) := sup |/i(/x) - t J log \T'\ dfjL-.fjLZ m\ , 

where A4 denotes the space of T— invariant probability measures and h(/j,) the 
entropy of the measure fx (see |Wa[ Chapter 4] for the definition and properties). 
For each tgR there exists a unique measure \it attaining the supremum in equation 
© that we call equilibrium measure (see |PU1 Chapter 5]). It is well known that 
the pressure function t — > P(t) is real analytic, strictly convex, strictly decreasing 
and that P'{t) = - JTog|T'| dfj, t (see [EE Chapter 3,4] and [HO Chapter 5]). 
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3.1. Newton map. A classical result by Bowen [Bo later generalised by Ruelle 
[Rul] establishes a relation between the topological pressure and the Hausdorff 
dimension of the repeller. 

Proposition 3.1 (Bowen's formula). Let T be a cookie cutter map and let d = 
dimn A then d is the unique root of the equation 

P(-tlog\T'\)=0. 

Let us denote by Np(t) the Newton map of the pressure function P(— t\og \T'\), 
that is 

P'(~t\og\T'\) 

We have that Np(d) = d and the iterations of Np(t) converge to the unique root 
of the pressure, that is 

lim Np-(t) = d. 

3.2. Lyapunov spectrum. The Lyapunov exponent of any given cookie cutter 
can attain an interval of values [a m - ln , a ma x] ■ For a £ [a m i n ,a max ] consider, 

J(a) = {x G / : \(x) = a} and J' = {x G / : the limit X(x) does not exists} . 

These level sets induce a decomposition of the repeller A of T. Indeed, A = 
J' U (U a J(a)) . The function that encodes this decomposition is called Lyapunov 
spectrum and it is defined by 

L(ct) := dim# ( J(a)). 

H. Weiss [We , building up on previous joint work with Pesin |PWj . proved that 
the Lyapunov spectrum is real analytic. A rather surprising result in light of the 
fact that this decomposition is fairly complicated. For instance, each level set turns 
out to be dense in A (see [We]). By means of the thermodynamic formalism it is 
possible to obtain formula for the Lyapunov spectrum. Indeed, 

(7) L(a) = iinf(P(-tlog|T'|)+ta). 

a teR 

This formula follows form the work of Weiss [We and can be found explicitly, for 
instance, in the work of Kessebohmer and Stratmann [KS (see also |IKj for some 
explicit computations). Actually, in this setting, the Lyapunov spectrum can be 
written as 

L{a) = -{P{-t a \og\T'\)+t a a), 
a 

where t a is the unique real number such that 

-±P(-t\og\T'\) 
at t=t a 

If we denote by [i a the unique equilibrium measure corresponding to the potential 
— t a log \T'\ we have that 



J log \ T '\ d^a = a. 



Therefore, 



(8) L(a) = -(P(-t a \og\T'\) + t a a), 

a 
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After the substitution a — a(t), equation (jHJl becomes: 

(9) L(a(t)) = (P(-t log |T'|) + ta(t)) . 



4. Main result 

In virtue of the comments we have made so far, we can now establish a relation 
between the Lyapunov spectrum of a cookie cutter map and the Newton map 
corresponding to the pressure function. 

Theorem 4.1. Let T be a cookie cutter map, then for every t S K we have 
(10) L(-P'(-t log \T'\)) = N P (t). 

Proof. The result follows noticing that the Lyapunov spectrum is given by (see 
equation ([9])) 

L(-P'(-t\og\T'\)) = _ p , { J loglT , l} (P(-tlog\T'\)-tP'(-tlog\T'\)). 
This combined with Lemma [5] yields the desired result. □ 

References 

[Ba] L. Barreira Dimension and recurrence in hyperbolic dynamics. Progress in Mathematics, 

272. Birkhuser Verlag, Basel, 2008. xiv+300 pp. 
[Bo] R. Bowen, Hausdorff dimension of quasi- circles. Inst. Hautes Etudes Sci. Publ. Math. 50, 

11-26 (1979). 

[Fa] K. Falconer Fractal geometry. Mathematical foundations and applications. Second edition. 
John Wiley and Sons, Inc., Hoboken, NJ, 2003. 

[IK] G. Iommi and J. Kiwi The Lyapunov spectrum is not always concave Journal of Statistical 
Physics vol.135 no.3, 535-546 (2009). 

[K] G. Keller, Equilibrium states in ergodic theory, London Mathematical Society Student Texts, 
42. Cambridge University Press, Cambridge, 1998. 

[KS] M. Kessebohmer and B. Stratmann A multifractal analysis for Stern-Brocot intervals, con- 
tinued fractions and Diophantine growth rates Journal fur die reine und angewandte Math- 
ematik (Crelles Journal) 605, 133-163, (2007). 

[Mi] J. Milnor Dynamics in one complex variable. Third edition. Annals of Mathematics Studies, 
160. Princeton University Press, Princeton, NJ, 2006. viii+304 pp. 

[PP] W. Parry and M. Pollicott, Zeta functions and the periodic orbit structure of hyperbolic 
dynamics. Astrisque No. 187-188 (1990), 268 pp. 

[PW] Y. Pesin, and H. Weiss A multifractal analysis of equilibrium measures for conformal ex- 
panding maps and Moran-like geometric constructions. J. Statist. Phys. 86 (1997), no. 1-2, 
233-275. 

[Pe] Y. Pesin Dimension Theory in Dynamical Systems CUP (1997). 

[PU] F. Przytycki and M. Urbahski, Conformal Fractals: Ergodic Theory Methods, Cambridge 
University Press 2010. 

[Rul] D. Ruelle Repellers for real analytic maps. Ergodic Theory and Dynamical Systems 2 
(1982), no. 1, 99107. 

[Ru2] D. Ruelle, Thermodynamic formalism. The mathematical structures of equilibrium statis- 
tical mechanics, Second edition. Cambridge Mathematical Library. Cambridge University 
Press, Cambridge, 2004. 

[Wa] P. Walters, An introduction to ergodic theory, Graduate Texts in Mathematics 79. Springer- 

Verlag, New York-Berlin, 1982. 
[We] H. Weiss The Lyapunov spectrum for conformal expanding maps and axiom-A surface dif- 

feomorphisms. J. Statist. Phys. 95 no. 3-4, , 615-632 (1999). 



6 



GODOFREDO IOMMI 



Facultad de Matematicas, Pontificia Universidad Catolica de Chile (PUC), Avenida 
Vicuna Mackenna 4860, Santiago, Chile 
E-mail address: giommiSmat .puc . cl 
URL: http: //www .mat .puc . cl/~giommi/ 



